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Abstract 

The Newton-Sabatier method for solving inverse scattering problem with fixed- 
energy phase shifts for a sperically symmetric potential is discussed. It is shown that 
this method is fundamentally wrong: in general it cannot be carried through, the 
basic ansatz of R.Newton is wrong: the transformation kernel does not have the form 
postulated in this ansatz, in general, the method is inconsistent, and some of the 
physical conclusions, e.g., existence of the transparent potentials, are not proved. A 
mathematically justified method for solving the three-dimensional inverse scattering 
problem with fixed-energy data is described. This method is developed by A.G. Ramm 
for exact data and for noisy discrete data, and error estimates for this method are 
obtained. Difficulties of the numerical implementation of the inversion method based 
on the Dirichlet-to-Neumann map are pointed out and compared with the difficulty 
of the implementation of the Ramm's inversion method. 

1 Introduction and conclusions 

It is shown in this paper that the Newton-Sabatier (NS) method for solving inverse scat- 
tering problem with fixed-energy phase shifts for spherically-symmetric potentials is fun- 
damentally wrong. The Ramm's method for solving inverse scattering problem with fixed- 
energy data is described for the exact data and for the noisy data, and the error estimates 
of this method are given. An inversion method using the Dirichlet-to-Neumann (DN) map 
is discussed, the difficulties of its numerical implementation are pointed out and compared 
with the difficulties of the implementation of the Ramm's inversion method. 

The inverse scattering problem (ISP) with fixed-energy data is of basic interest in 
quantum and classical physics and in many applications. The statement of the problem 
is well known and can be found in [8]. In [4] and [2] the NS method for solving ISP for 
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spherically symmetric potentials is described. In the sixties P. Sabatier published several 
papers concerning this procedure, and there are quite a few papers of several authors using 
this procedure and generalizing it. A vast bibliography of this topic is given in [2] and [4], 
and by this reason we do not include references to many papers treating this topic. 
The NS method for finding q{r),r := \x\, consists of the following [2], [4]-[6]: 

{Se}i=o,i,2, {ci} Kir, s) ^ g^(r) := — ^ (1-1) 

Here 6e are the fixed-energy phase shifts, generated by some potential q{r), q are some 
constants which should be calculated from 6^ by solving an infinite linear algebraic system, 
the K{r, s) has to be found from the equation 

K{r, s) = fir, s) - f Kir, s)*-^ dt, (1.2) 

JO 

where 

fir,s) ■.= 22(^eMr)Ms), utikr) := -J —Jt+i/2ikr), (1.3) 

and is the Bessel function. 

Our conclusion is: NS method is fundamentally wrong and, in general, 
cannot be carried through. In the exceptional cases, when it can be carried 
through, it does not yield the original generic potential which generated the 
data {(5£}f=o,i,2, • 

By a generic potential q we mean a q, which is not a restriction to (0, oo) of an analytic 
function. 

It is not proved in [4] and [2] that the qj^ from (1.1) generates the original data 
{^^}£=o,i,2, • 

The R.Newton's ansatz (1.2)-(1.3) for the transformation kernel K(r, s) is wrong: in 
general, Kir,s) does not solve (1.2)-(1.3). 

The set of potentials q^ with YIILq |q| < oo, is not dense in Li^i '■= {q q = q, \\q\\ < 

Ikll •= /o°° ''kl'^^}' ^ standard scattering class. Therefore the NS method cannot be 
used even as a parameter-fitting procedure for solving ISP. 

A detailed discussion and a justification of the above conclusions is given in [9] and 
[12] (see also [1]). For convenience of the reader we give a brief justification here. 

First, R.Newton did not prove that eq. (1.2) is solvable for all r > 0. If it is not 
solvable for at least one r > 0, then the NS method breaks down: it yields q^ which is not 
locally integrable [9], and the scattering theory cannot be constructed for such potentials, 
in general. 

Secondly, R. Newton did not prove existence of the transformation kernel Kir,s), 
independent of £ and did not study its properties. This was done in [10] , and it turns out 

that R. Newton's ansatz (1.2)-(1.3) for K, which was not justified by R. Newton in any 
way, is wrong. In particular, if K solves (1.2)-(1.3), and q are growing not faster than 
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an exponential, then the function K{r, r) must be a restriction to (0, oo) of an analytic 
function, which is not the case for a generic potential. 

In the above I assumed that the inverse scattering problem, which R. Newton tried to 
solve, consists of finding a q belonging to a certain functional class, (for example, to Li^i, 
a standard scattering class), given the fixed-energy phase shifts generated by this q, and 
I have demonstrated that the NS method does not solve this problem. 

In fact, the NS method does not solve even a much less interesting problem: given the 

data {5^}^=o,i,2, ) corresponding to some potential q, find a potential qn which generates 

the same data. This problem is much less interesting than the inverse scattering problem, 
because q^ may have no physical meaning, for example, it may deacay very slowly, etc. 
The NS method does not solve even this problem because equation (1.2) may be not 
solvable for some r > 0, and also because there is no proof that gjv generates the original 
data {de}e=o,i,'2, • 

In [2] the existence of the "transparent" potentials is claimed, that is, potentials which 
produce, at a fixed energy, all the phase shifts equal to zero. This claim is not proved 
because there is no proof of the existence of the solution of the corresponding equation 
(1.2) for all r > 0, and also because there is no proof that q^ generates the original data 
{Se}e=o,i,2, ■ 

Let us finally explain why the set of potentials qjsf, which can possibly be obtained by 
the scheme (1.1) with l*-^! ^ dense in Li^i in the norm || • ||. Note that 

in the examples, considered in [4] and [5], q = for all sufficiently large I, so that the 
condition YltLo |q| < oo is satisfied. 

Assume that q G Li^i and that there is a sequence of q^j, denoted qj to simplify the 
notation, such that — ^ as j ^ cxd. Define L{q) := rqdr. Clearly L is a linear 
functional on Li^i and ||L|| = 1. Thus L{qj) L{q) as j oo. Choose q with L{q) ^ 0. 
Then L{qj) ^ for sufficiently large j. One has -fir(r, r) = — | sqN{s)ds, and therefore 



because in our case L{qj) L{q) ^ 0. Thus, one has a contradiction, because [^(r, r)| < 
sup^>oS^o k^ll^^('")llw('')l < where ipi{r) is the regular solution of the Schroedinger 
equation ip'^ + ip£ — ^^^V^ — Q{f)ip^ = 0, r > 0, </? ~ (2^+1)!! as r — >■ 0. Here we have 
used the relations: 



where 5^ are the phase shifts at = 1, and is the Jost function at = 1. It can be 
proved that sup^ \F(\ < 00. Thus, if YleLo k^l < then K(r,r) = 0{1) as r — ^ 00. 
In [12] a detailed analysis of the NS method is given. 

In [3] a uniqueness theorem is claimed for an equation which is a version of (1.2). 
In [11] a counterexample to this wrong theorem is constructed. In Section 2 Ramm's 
inversion method is described. In Section 3 an inversion method which uses the DN map 
is described. 



K{r, r) = — 



[1 + 0(1)] — 00 as r — 00, 
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2 Ramm's inversion method for exact data 



The results wc describe in this Section are taken from [8] and [13]. Assume q G Q := 
Qa n L°°(M3), where Qa := {q : = ^(.t) G L^{Ba), q{x) = if \x\ > a}, 

Ba := {x : \x\ < a}. Let A(a',a) be the corresponding scattering amphtude at a fixed 
energy A;^, A; = 1 is taken without loss of generality. One has: 

oo „ 

A{a', a)=J2 Ae{a)Ye{a'), Ae{a) := / ^ A{a', a)Yei^da' , (2.1) 

where 5^ is the unit sphere in R^, Y£{a') = Y£^m{a'), —£ < m < £, are the normalized 
spherical harmonics, summation over m is understood in (2.1) and in (2.8) below. Define 
the following algebraic variety: 

3 

M ■-{e:9GC^,e-e = i}, e -w — ^ejWj. (2.2) 

This variety is non-compact, intersects over S'^, and, given any ^ G R^, there exist 
(many) 0,0' e M such that 

0'-0 = ^^ 1^1 ^oo, 0,0' eM. (2.3) 

In particular, if one chooses the coordinate system in which ^ = te^, t > 0, 63 is the unit 
vector along the X3-axis, then the vectors 

= ^e3 + C2e2 + Ciei, ^ = -^eg + C2e2 + Ciei, C? + Cl = l-j, (2.4) 

satisfy (2.3) for any complex numbers and C2 satisfying the last equation (2.4) and such 
that ICiP + IC2P — OO- There are infinitely many such (1,(2 € C. Consider a subset 
M' C M consisting of the vectors 9 = (sini? cos (^9, sini? sin 99, cos 1? where 'd and ip run 
through the whole complex plane. Clearly G M, but M' is a proper subset of M. Indeed, 
any e M with 0^ 7^ ±1 is an element of M'. If ^3 = ±1, then cosi? = ±1, so sini? = 
and one gets = (0, 0, ±1) G M'. However, there are vectors = {0i,02,l) € M which 
do not belong to M' . Such vectors one obtains choosing ^i, ^2 S C such that 9\ + 02 = 0. 
There are infinitely many such vectors. The same is true for vectors {0i,02,—^)- Note 
that in (2.3) one can replace M by M' for any ^^M?, ^^2e^. 
Let us state two estimates proved in [8]: 

^\AAa)\<.{^y{^f\ (2.5) 

and 

1 pr\Imd\ 

|y,(0)| < -=-^, Vr>0, 0eM', (2.6) 
V47r \uir]\ 
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where 

Mr) := (^) ' (r) = ^] (|)' [1 + o(l)] as £ ^ oo, (2.7) 

and Je{r) is the Bessel function regular at r = 0. Note that Y£{a'), defined above, admits 
a natural analytic continuation from S'^ to M by taking ?? and if to be arbitrary complex 
numbers. The resulting 9' € M' C M. 

The series (2.1) converges absolutely and uniformly on the sets S"^ x Mc, where Mc is 
any compact subset of M. 

Fix any numbers oi and b, such that a < ai < b. Let || • || denote the L^{ai < \x\ < b)- 
norm. If \x\ > a, then the scattering solution is given analytically: 



uix,a) =e'"'' + ^Ai{a)Ye{a')heir), r := \x\ > a, «':=-, (2.8) 

where A£{a) and y^(a') are defined above, 

heir) :=e^f(^+i)./^ff« (r), 

i?^^^(r) is the Hankel function, and the normalizing factor is chosen so that /i^(r) = 
^[1 + o(l)] as r —>■ oo. Define 

p{x) := p{x; u) := e"^^'"^ / ^(a;, a)z/(a, e)da - 1, G l2(52). (2.9) 
Consider the minimization problem 

IIpII = inf := d{e), (2.10) 

where the infimum is taken over all v G L^(S'^), and (2.3) holds. 
It is proved in [8] that 

d{e) < c\e\-^ ifOeM, \e\ » i. (2.11) 

The symbol |^| » 1 means that |^| is sufficiently large. The constant c > in (2.11) 
depends on the norm ||q'||L2(Ba) but not on the potential q{x) itself. 

An algorithm for computing a function i/(a, 6), which can be used for inversion of the 
exact, fixed-energy, three-dimensional scattering data, is as follows: 

a) Find an approximate solution to (2.10) in the sense 

\\p{x,,^)\\<2d{9), (2.12) 
where in place of the factor 2 in (2.12) one could put any fixed constant greater than 1. 
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b) Any such z^(a, 6*) generates an estimate of q{$,) with the error O ^]^^) 1^1 ~^ 
This estimate is calculated by the formula 

q:=-4.Tr [ A{e',a)u{a,e)da, (2.13) 

where v{a,6) G L?'{S'^) is any function satisfying (2.12). 
Our basic result is: 

Theorem 2.1. Let (2.3) and (2.12) hold. Then 

sup|g-^(6l<TZT' I^Hoo, (2.14) 

The constant c > in (2.14) depends on the norm of q, but not on a particular q. 

In [8] and [13] an inversion algorithm is formulated also for noisy data, and the error 

estimate for this algorithm is obtained. Let us describe these results. 

Assume that the scattering data are given with some error: a function As {a', a) is 
given such that 

sup \A{a',a)-As{a',a)\<6. (2.15) 

a',ae52 

We emphasize that As{a' , a) is not necessarily a scattering amplitude corresponding to 
some potential, it is an arbitrary function in L°°{S'^ x S*^) satisfying (2.15). It is assumed 
that the unknown function A{q! , a) is the scattering amplitude corresponding to a g G Q. 

The problem is: Find an algorithm for calculating qs such that 

sup \qs-q{0\<riiS), r?(<5) ^ as 5 ^ 0, (2.16) 

and estimate the rate at which r]{S) tends to zero. 

An algorithm for inversion of noisy data will now be described. 
Let 

llnJI 



N{6): 

where [x] is the integer nearest to x > 0, 



ln|ln(5|J ' 



(2.17) 



N{5) 

Asie',a):=Y,^Siia)Yi{e'), Asi{a) := As{a' ,a)Yt{a')da' , (2.18) 



N{5) 



us{x, a) := e^"-^ + ^ Ase{a)Ye{a')he{r), (2.19) 

ps{x; v) := e-'^"" / us{x, a)v{a)da - 1, 6* G M, (2.20) 
^52 
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^(S) := e-T^(^) , 7 = In ^ > 0, (2.21) 

a(zy) := ||i/||^2(52), k := \Ime\. (2.22) 
Consider the variational problem with constraints: 

1^1 = sup := T^{d), (2.23) 

1^1 [\\psii^)\\ + a(i.)e'^V('^)] < c, 9eM, \e\ = sup := ^{S), (2.24) 

the norm is defined above (2.8), and it is assumed that (2.3) holds, where ^ G is an 
arbitrary fixed vector, c > is a sufficiently large constant, and the supremum is taken 
over 6 e M and u G L?{S'^) under the constraint (2.24). By c we denote various positive 
constants. 

Given ^ G one can always find and 6' such that (2.3) holds. We prove that 
'd{5) — DO, more precisely: 

Let the pair 6{5) and Uj{a,6) be any approximate solution to problem (2.23)-(2.24) in 
the sense that 

m 



m\ > (2.26) 



Calculate 



qs-'-A-K As{e',a)i^sia,9)da. (2.27) 

Theorem 2.2. // (2.3) and (2.26) hold, then 

^, (ln|ln(5|)2 
sup|96-g(OI <c I .1 as5^Q, (2.28) 
^gjjs |ind| 

where c > a constant depending on the norm of q. 

In [8] estimates (2.14) and (2.28) were formulated with the supremum taken over an 
arbitrary large but fixed ball of radius ^o- Here these estimates are improved: ^0 = 00. 
The key point is: the constant c > in the estimate (2.11) does not depend on ^. 

Remark. In [14] (see also [15] and [13]) an analysis of the approach to ISP, based on 
the recovery of the DN (Dirichle-to-Neumann) map from the fixed-energy scattering data, 
is given. This approach is discussed in Section 3. 

The basic numerical difficulty of the approach described in Theorems 2.1 and 2.2 
comes from solving problems (2.10) for exact data, and problem (2.23)-(2.24) for noisy 
data. Solving (2.10) amounts to finding a global minimizer of a quadratic form of the 
variables if one takes u in (2.9) as a linear combination of the spherical harmonics: 
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V = X]^=o '^^^('^)- If o'^^ uses the necessary condition for a minimizer of a quadratic form, 
that is, a Unear system, then the matrix of this system is ih-conditioned for large L. This 
causes the main difficulty in the numerical solution of (2.10). On the other hand, there 
are methods for global minimization of the quadratic functionals, based on the gradient 
descent, which may be more efficient than using the above necessary condition. 



3 Discussion of the inversion method which uses the DN 
map 

In [14] the following inversion method is discussed: 

q{i)= lim e^p{-ie' ■ s){A- Ao)i;ds, (3.1) 

where (2.3) is assumed, A is the Dirichlet-to-Neumann (DN) map, ijj is found from the 
equation: 

V'(s) = ipois) - G{s- t)Bipdt, B :=A- Aq. (3.2) 
Js 

and G is defined by the formula: 

The DN map is constructed from the fixed-energy scattering data A(a', a) by the method 
of [14] (see also [8]). 

Namely, given A{a',a) for all a', a G S'^, one finds A using the following steps. 

Let / G H^/^{S) be given, S' is a sphere of radius a centered at the origin, are its 
Fourier coefficients in the basis of the spherical harmonics, 

w = f2fiYi{x^)^,, r>a, :=-, r:=\x\. (3.4) 



Let 

w 



/ g{x,s)a{s)ds, (3.5) 
Js 



where a is some function, which we find below, and g is the Green function (resolvent 
kernel) of the Schroedinger operator, satisfying the radiation condition at infinity. Then 

wIj = w^ + cr, (3.6) 

where N is the outer normal to 5", so N is directed along the radius- vector. We require 
w = f on S. Then w is given by (3.4) in the exterior of S, and 

•"^ = |/.^.(x»)^. (3.7) 
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By formulas (3.6) and (3.7), finding A is equivalent to finding a. By (3.5), asymptotics of 

w as r := \x\ oo, x/\x\ := x^, is (cf [8, p. 67]): 

e''' u{y,-x^) 1 

w = h o(-), (3.8) 

r 47r ^ ' 

where u is the scattering solution, 

oo 

u{y, -x^) = e— + ^ Ae{-x^)Ye{y^)he{\y\). (3.9) 
e=o 

From (3.4), (3.8) and (3.9) one gets an equation for finding a ([14, eq. (23)], see also [8, 
p. 199]): 

which can be written as a linear system: 

-y^ = a2(-l)' ^aH47ri'jK«)'^"zF + A'//i«'(a)], (3.11) 

''^^^ l'=0 

for the Fourier coefficients of o". The coefficients 

are the Fourier coefficients of the scattering amplitude. Problems (3.10) and (3.11) are 
very ill-posed (sec [14] for details). 

This approach faces heavy difficulties: 

1) The construction of the DN map from the scattering data is a very ill-posed problem, 

2) The construction of the potential from the DN map is a very difficult problem 
numerically, because one has to solve a Fredholm-type integral equation ( equation (3.2) 
) whose kernel contains G, defined in (3.3). This G is a tempered distribution, and it is 
very difficult to compute it, 

3) One has to calculate a limit of an integral whose integrand grows exponentially to 
infinity if a factor in the integrand is not known exactly. The solution of equation (3.2) 
is one of the factors in the integrand. It cannot be known exactly in practice because it 
cannot be calculated with arbitrary accuracy even if the scattering data are known exactly. 
Therefore the limit in formula (3.1) cannot be calculated accurately. 

No error estimates are obtained for this approach. 

In contrast, in Ramm's method, there is no need to compute G, to solve equation 
(3.2), to calculate the DN map from the scattering data, and to compute the limit (3.1). 
The basic difficulty in Ramm's inversion method for exact data is to minimize a quadratic 
form (2.10), and for noisy data to solve optimization problem (2.23)-(2.24). The error 
estimates are obtained for the Ramm's method. 
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